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We report the full complex dielectric function of high-purity Bao.68Ko.32Fe2As2 single crystals with Tc — 
38.5 K determined by wide-band spectroscopic ellipsometry at temperatures 10 < T < 300 K. We discuss 
the microscopic origin of superconductivity-induced infrared optical anomalies in the framework of a multiband 
Eliashberg theory with two distinct superconducting gap energies 2Aa ~ 6 fcaTc and 2Ab ~ 2.2 /cbTc. The 
observed unusual suppression of the optical conductivity in the superconducting state at energies up to 14 fcaTc 
can be ascribed to spin-fluctuation-assisted processes in the clean limit of the strong-coupling regime. 
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The discovery of iron-based superconductors [1] has gen- 
erated significant experimental and theoretical effort to un- 
ravel the mechanism of high-temperature superconductivity 
in these compounds. This effort has yielded a comprehen- 
sive experimental description of the electronic structure at the 
Fermi level, which includes multiple Fermi surface sheets in 
a good agreement with density functional calculations [2, 3]. 
Partial nesting between at least two of these sheets leads to a 
spin-density-wave instability that renders the metallic parent 
compounds antiferromagnetic. In the superconducting com- 
pounds spin fluctuations become the source of strong repul- 
sive interband interactions and might give rise to supercon- 
ductivity with different signs on these sheets [4]. 

The most incisive experimental data have been obtained 
on high-quality single-crystals of iron pnictides with the so- 
called 122 structure, for instance BaFe2As2, with K substi- 
tuted for Ba or Co for Fe, resulting in hole and electron dop- 
ing, respectively. In all of these materials five Fermi surface 
sheets have been identified in calculations and confirmed by 
numerous independent experimental studies [2, 3]: in the re- 
duced Brillouin-zone scheme these are three hole pockets at 
the r point and two almost degenerate electron pockets at 
the X point with nesting between hole and electron sheets. 
Among all of these 122 materials, the optimally hole-doped 
compound Bao.68Ko.32Fe2As2 (BKFA) has the highest tran- 
sition temperature of 38.5 K. Due to their exceptional quality, 
crystals of this compound are well suited as a testbed for the- 
oretical models. A four-band Eliashberg theory with strong 
interband couping has akeady proven successful in account- 
ing for the transition temperature, as well as the temperature 
dependence of the free energy and superconducting gaps of 
this compound [5]. This analysis has made clear that a satis- 
factory description of the bulk thermodynamical properties in 
the superconducitng state can only be obtained via strong cou- 
pling to spin fluctuations or other bosons with spectral weight 
below 50 meV. 

In this Letter we extend the approach of Ref. 5 to describe 
the far-infrared properties of the same BKFA single crystal. 
We show that major characteristic features of superconduc- 
tivity can be explained within a strong-coupling Eliashberg 
approach with two distinct values of the superconducting en- 
ergy gap 2Aa ~ 6 k^Tc and 2Ab ~ 2.2 k^T^, in quan- 
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FIG. 1: Real part of the (a) optical conductivity and (b) dielectric 
function in the far-infrared spectral region. Two characteristic super- 
conductivity energy scales are present: 6 fcaTc and 14 fceTc. 



titative agreement with angle-resolved photoemission [6-8], 
scanning-tunneling microscopy [9] and specific-heat mea- 
surements [5]. We also demonstrate that within this approach 
the qualitative differences in the infrared spectra of electron- 
and hole-doped 122 compounds are reproduced by strong- 
coupling calculations in the clean and dirty limits (weak and 
strong impurity scattering), respectively. 

The optimally-doped BKFA single crystals were grown 
in zirconia crucibles sealed in quartz ampoules under argon 
atmosphere [10]. From DC resistivity, magnetization and 
specific-heat measurements we obtained Tc = 38.5 ± 0.2 K. 
The sample surface was cleaved prior to every optical mea- 
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FIG. 2: Optical scattering rate obtained from the experimental data 
at 41 and 10 K within the extended-Drude model, with the contri- 
bution of the interband transitions subtracted (blue and black lines, 
respectively) and at 41 K without subtraction (green line). Dash- 
dotted line indicates the saturation level of the high-energy optical 
scattering rate. 



surement. The full complex dielectric function e{uj) was ob- 
tained in the range 0.01 — 6.5 eV using broadband ellipsome- 
try, as described in Ref. [11]. In this work we focus on the itin- 
erant charge carrier response contained within the far-infrared 
spectral range measured at the infrared beamline of the ANKA 
synchrotron light source at Karlsruhe Institute of Technology, 
Germany. The contribution of the interband transitions has 
been eliminated based on a dispersion analysis in the entire 
spectral range [12]. 

The optical response of BKFA in the far-infrared spectral 
range is shown in Fig. l(a, b) respectively for the real parts 
of optical conductivity (T(aj) — cri{Ld) + i<T2{Ld) and dielec- 
tric function e{oj) — 1 + 47ri(T(tj)/cj. It is dominated by 
the contribution of the itinerant charge carriers manifested 
in negative values of £i(w). Figure 1(a) also reveals a peak 
around 15 — 20 meV in cri(w) with a concomitant upturn 
in ei{uj) at higher temperatures indicating the presence of a 
collective excitation. The superconducting transition is ac- 
companied by the suppression of the optical conductivity up 
to 50 meV (14fcBTc). The corresponding missing area in 
cri(w) between 41 and 10 K, Acri(w)da; = (8A£)-i, 

manifests itself as a characteristic — 1/(Alw)^ contribution to 
ei{uj) in Fig. 1(b) in the superconducting state. The London 
penetration depth Al = 2200 A extracted from these data is 
consistent with other measurements [13]. At energies close to 
the optical superconducting gap 2Aa ~ 20 meV one of the 
directly measured ellipsometric angles \l/(a;) approaches its 
critical value of 45° at the superconducting transition, which 
implies that the reflectivity of the sample approaches unity and 
its optical conductivity cri(w) is close to zero. Remarkably, 
Fig. 1(a) shows a quasilinear dependence of ai{uj) in the su- 
perconducting state from 2Aa to as high as lAk^Tc, in a stark 
contrast to the electron-doped 122 compounds [14-16]. In the 
latter the optical conductivity at 2Aa decreases abruptly upon 
cooling below Tc, but only a weak superconductivity-induced 



modification is observed at higher energies. The quasilinear 
behavior in BKFA cannot be reconciled with the widely used 
for pnictides Mattis-Bardeen theory [17], a weak-coupling ex- 
tension of the BCS theory to finite impurity scattering. As all 
optimally-doped 122 pnictide superconductors appear to be in 
the strong-coupling regime, the Eliashberg theory [18] has to 
be used in order to obtain an adequate description of the opti- 
cal properties. 

Signatures of a boson pairing mediator of the Eliashberg 
theory come from a qualitative analysis of the optical con- 
ductivity within the extended Drude model. It implies that 
the optical scattering rate is related to the far-infrared opti- 
cal response as ^{oj) = Re[a;pj/47r(7'*(a;)], where the su- 
perscript 'it' stands for 'itinerant' and implies that the con- 
tribution of all interband transitions must be subtracted from 
the experimentally obtained optical conductivity. The optical 
scattering rate shows clear evidence of an intermediate bo- 
son irrespective of complications due to the multiband char- 
acter of the compound. Figure 2 plots 7(0;) of BKFA at 
41 K (blue line) and 10 K (black line) for Wpi = [SSWit]^/^ ^ 
[8 /p°°cri'(w)dw]i/2 = 1.6 eV with afl interband transitions 
subtracted in both cases (SWit corresponds to the blue shaded 
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FIG. 3: (a) Real part of the optical conductivity at 41K (blue line). 
The contribution of itinerant charge carriers (blue area) is obtained 
by subtracting all interband transitions ai"''''^(cj) (gray area) from the 
optical response, (b) Real part of the dielectric function at 41K (blue 
line). The free-charge-carrier response ei'(i^) (open circles) is ob- 
tained by eliminating all interband transitions e^"*'^'^ [uj) (black solid 
line). The blue dashed line indicates the screened plasma frequency 
at41K. 
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area in Fig. 3(a)). In the superconducting state, no scattering is 
expected up to photon energies exceeding the binding energy 
of the Cooper pairs. Thus the onset of the optical scattering 
rate marks the optical energy gap 2Aa = 20 meV. Satura- 
tion of 7(w > 50 meV) at 1100 cm^^ indicates that the boson 
spectral function is contained well below 50 meV [19]. It is 
important to emphasize that, due to the multiband character 
of the iron pnictides, an analysis of the optical scattering rate 
in the framework of a single-band Eliashberg theory is poten- 
tially misleading. Moreover, also shown in Fig. 2 is a spec- 
trum that directly results from the experimental data, without 
accounting for the interband transitions. It becomes clear that 
an increase in the scattering rate at higher energies that might 
be ascribed to strong electron correlations can result from an 
unsubtracted contribution of the interband transitions to the 
complex optical conductivity. This is especially important in 
iron pnictides since the lowest lying interband transition at 
about 0.5 eV contributes to an anomalously large value of the 
low-energy dielectric permittivity Eoo [20] due to the high po- 
larizability of the Fe-As bonds [21]. In order to reconcile the 
bare plasma frequency of 1.6 eV (see Fig. 3(a)) with the zero- 
crossing in £1(0;) at 0.2 eV (blue line in Fig. 3(b)) £^0 has to 
be as large as 60, consistent with the contribution of the in- 
terband transitions e"*''' (a;) determined by means of the dis- 
persion analysis, as shown in Fig. 3(b). Such Soo is thus an 
order of magnitude larger than in any other high-temperature 
superconductor (e. g. « 5 in cuprates [22]). Recently, a simi- 
larly high value in a conventional superconductor was inferred 
from reflectivity measurements on elementary bismuth [23]. 

To determine the microscopic origin of the high-energy 
anomaly 2Aa < hu < 14fcBrc in the real part of the op- 
tical conductivity in Fig. 1(a) we use a four-band Eliashberg 
theory that proved successful in explaining thermodynamical 
data obtained on the same compound [5]. The 4x4 matrix 
of coupling constants and 4 densities of states characterizing 
this model are highly constrained by thermodynamic, trans- 
port and photoemission data [5-9], and the same set of param- 
eters is used here. In principle, an additional set of 4 plasma 
frequencies and a 4 x 4 matrix of intraband/interband impu- 
rity scattering rates has to be taken into account to describe the 
optical response. However, this parameter set can be strongly 
reduced based on the following considerations. 

A substantial simplification is made possible by a projec- 
tion of the four-band model onto an effective two band model 
motivated by the observation of two distinct groups of su- 
perconducting energy gaps in a variety of experiments [5- 
9]. These gaps can be identified as a single gap As on the 
outer hole-like Fermi surface and a group of three gaps of 
magnitude ~ A^ on the inner hole-like and the two electron- 
like Fermi surfaces. Minimizing the ground-state energy sub- 
ject to this grouping constraint yields an effective two-band 
model [12]. Furthermore, as the superconducting transition 
temperature of BKFA appears to be only weakly correlated 
with the residual resistivity (which is a measure of the impu- 
rity scattering), off-diagonal elements of the impurity scatter- 
ing matrix can be neglected (see Table SI in Ref. 12). 




20 40 60 

Photon energy (meV) 



FIG. 4: (a) Real part of the far-infrared conductivity obtained within 
the two-band Eliashberg theory (see text) at 40 K (blue lines) and 
10 K (black lines) in the clean limit 7a = 7b = 1 cm~^ (solid lines) 
and dirty limit 7a = 7b = 200 cm~^ (dashed lines), (b) Optical 
scattering rate in the clean limit from the same model. The gray 
area shows the normalized boson spectral function B{ij) used in the 
calculation, displaced from zero by 2Aa to assist interpretation in 
the superconducting state [25]. 



Given the boson spectrum centered at 13 meV (see supple- 
mentary online material in Ref. 5) consistent with the energy 
of the spin resonance excitation in this compound [24, 25] 
one obtains the following two-band model coupling matrix: 
Aaa = 4.36, Abb = 0.2, Aab = -0.35, Aba = -0.5, 
with the fractional density of states being Na/{Na + Nb) = 
0.59 [12]. The first effective intraband coupling constant is 
an order of magnitude larger than predicted for the intraband 
electron-phonon coupling [26]. It does not, however, bear any 
physical meaning by itself but rather incorporates contribu- 
tions from three different bands. We reiterate that the cou- 
pling matrix has been inferred from prior measurements. In 
this way, only two intraband impurity scattering rates enter as 
free parameters of the theory in addition to the plasma fre- 
quencies of the bands. 

In our calculation we consider two clean bands with 7a = 
7b = 1 cm^^. As the bare plasma frequencies of all bands 
are similar it follows that the spectral weight of band A has 
to be much larger than that of band B. Assigning 80% of the 
spectral weight to the effective band we obtain the results pre- 
sented as solid lines in Fig. 4(a). The high-energy anomaly at 
14fcB7c is naturally captured by the model without resorting 
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to additional gaps. Its energy is given by 2Aa + where 
il is the characteristic frequency of the boson spectrum, as 
shown in Fig. 4(b) (gray shaded area, displaced from zero by 
2Aa). This calculation also accounts for the fact that only 
the biggest superconducting gap is visible in the optical re- 
sponce of BKFA due to a small contribution of band B (20% 
of the spectral weight) to the overall optical conductivity. This 
leads to two possible levels of the impurity scattering rate of 
band B, which has to be either very small 73 w 1 cm^^ or 
very large at about 1000 cm^^. The latter value provides a 
better description of the optical scattering rate (see interac- 
tive simulation in [12]) and DC transport [27]. However, such 
a large disparity between the charge carriers is hard to rec- 
oncile with the Hall and de Haas-van Alphen experiments, 
which imply that the impurity scattering rate of the holes is 
no more than one order of magnitude higher than that of the 
electrons [28, 29]. This residual uncertainty notwithstanding, 
our results show that the impurity scattering rate of band A 
must be very small, because the energy 2Aa + is no longer 
discernible in the simulated spectra when 7a increases (inter- 
active simulation in [12]). The region of linear increase of 
cri(aj) is related to the linear segment in the boson spectrum 
and can only be observed in a very clean material. 

The same reduced two-band model can be applied to the 
case of BaFei.85Coo.i5As2 (BFCA). In this compound the 
spin resonance excitation occurs at a very similar energy of 
10 meV [30]. A boson spectrum centered at this energy is 
also consistent with Andreev-reflection measurements [31]. 
Recently, a comprehensive specific-heat study of this com- 
pound at different Co-doping levels has been carried out [32]. 
The analysis of the experimental data in the framework of 
the two-band a-model indicates that the largest gap develops 
in the band with the largest electronic density of states, pro- 
viding further evidence that several bands contribute to the 
strongly-coupled band in the reduced two-band model. Fig- 
ure 4(a) (dashed lines) shows that a calculation within the 
same reduced two-band model qualitatively reproduces the 
far-infrared optical conductivity of BFCA [14-16] when both 
bands are assumed to be dirty with 7a = 7b = 200 cm^^ and 
a redistribution of the spectral weight between the bands is 
taken into account as ^ « 3. The model captures the 
two prominent superconductivity-induced anomalies clearly 
observed in experiments: the steep onset of absorption at the 
value of the small gap 2Ab and the weaker superconductivity- 
induced changes of the optical conductivity extending up to 
ISfceTc. The redistribution of the spectral weight between 
the bands in BFCA compared to BKFA implied by our anal- 
ysis is justified by doping with different carriers in the two 
compounds, while the large difference in the their impurity 
scattering rates is a natural consequence of the difference in 
doping mechanisms by chemical substitution, which directly 
affects the FeAs layers in BFCA, but not in BKFA. 

In summary, a qualitative description of superconductivity- 
induced optical anomalies in the far-infrared optical conduc- 
tivity of Bao.68Ko.32Fe2As2 is obtained in the framework of 
an effective two-band Eliashberg theory with a strong cou- 



pling to spin fluctuations reduced from its four-band counter- 
part. The linear increase of absorption above the larger su- 
perconducting gap can only be observed when the effective 
band is extremely clean. The same model in the dirty limit 
provides a good qualitative explanation of the optical conduc- 
tivity of the optimally electron-doped BFCA consistently in 
the strong-coupling regime. 
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DISPERSION ANALYSIS 

The full complex dielectric function e{uj) — e\ (w) +i£2 iy) 
obtained experimentally in the range from 12 meV to 6.7 eV 
was analyzed in the Drude-Lorentz model: 



£(C.) = 1 - 



pi 



where (wpi, 7) are the plasma frequency renormalized by in- 
teraction with the mediating boson and optical scattering rate, 
and (Asj, woj, Fj) are the DC permittivity contribution, cen- 
ter frequency and the width of the Lorentzian oscillators used 
to model the interband transitions, respectively. The results 
of this analysis are presented in Fig. SI for 10 K (blue line 
- experimental data, black lines - separate Lorentz contribu- 
tions). To display the scale of the temperature-induced varia- 
tion of the interband transitions the experimental spectrum at 
300 K is also shown (red line). The lowest interband transi- 
tion in this material lies around 0.5 eV and significantly con- 
tributes to the AC polarizability of the system, as is evident 
from Fig. 3(b) of the main text. The residual optical response 
(open circles and inset in Fig. 3(b)) was studied after the sub- 
traction of all thus determined interband transitions down to 
0.5 eV. An unscreened bare plasma frequency of 1.6 eV at 
41 K was consistently obtained from the spectral weight of 
the residual response SW = ^'^ af{uj)duj as ojp\ ~ \/8SW 
and a simultaneous fit of the real and imaginary parts of the 
dielectric function at high energies. 



EFFECTIVE TWO-BAND MODEL FOR PNICTIDES 
Multiband Eliashberg formalism 

The Eliashberg theory of superconductivity [S1,S2] ex- 
tended to the multiband case [S3] has already proven suc- 
cessful in describing the thermodynamical properties of iron 
pnictides [S4]. Here we employed the same formalism to 
account for the far-infrared optical response of these com- 
pounds. The main ingredient of the theory is the total spectral 
function of the electron-boson interaction B{uj) (Eliashberg 
function; analogous to that of the electron-phonon interaction 
a^F(w)). In a four band system it can be decomposed into 16 
functions B{uj)ij, where i and j label the four Fermi surface 
sheets (i,j — 1,2,3,4). The sfawt/arc/ Eliashberg functions 
determine the superconducting and thermodynamical proper- 
ties such as the superconducting transition temperature and 
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FIG. SI: Real part of the (a) optical conductivity and (b) dielectric 
function (b) at 300K (red line) and lOK (blue line). Interband transi- 
tions inferred from the dispersion analysis (gray lines). The shaded 
region is enlarged in Fig. 3(b) of the main text. 



gaps, electronic specifc heat, de Haas-van Alphen mass renor- 
malizations etc. and are defined as 



■'k-k 



k,k',iy 



where Ni is the partial density of states per spin on the i'th 
sheet of the Fermi surface, and g^^^, is the matrix element of 
electron-boson interactions. 

Transport and electrodynamical properties are defined by 
16 transport Eliashberg functions (which enter the Boltzmann 
kinetic equation) 



1 



2N,{v 



a 2\ 
Fi I 



E 



5k,k' 



k.k' 



where is the a-th Cartesian component of the Fermi veloc- 
ity on Fermi surface i. The average Fermi velocity is related 
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TABLE S 1 : Superconducting transition temperature and residual re- 
sistivity P40K of (nearly) optimally hole-doped BKFA. 



Tc, K Residual resistivity, mQ cm 


Reference 


38.5 0.04 


[S4] 


38 0.075 


[S7] 


38 0.1 


[S6] 


36.5 0.055 


[S8] 



to the plasma frequency by the standard expression w^j ^ = 
STre^iY, (vl,) = Sttg^ vl,{k)S{el). All Eliashberg func- 
tions satisfy the symmetry relations MiBij — MjBji, where 
Mi — Ni and = w^jj for the standard and transport 
EUashberg functions, respectively. 



Role of impurities and defects 

Both normal and superconducting properties of a multiband 
superconductor significantly depend on impurity scattering. 
Unlike in conventional superconductivity, one has to distin- 
guish between the intraband impurity scattering, which does 
not add any new physics (in the Born approximation) com- 
pared with single-band superconductivity, and interband scat- 
tering, which in many cases has an effect comparable to the 
pair-breaking effect of magnetic impurities (or of nonmag- 
netic impurities in superconductors with p- or d-wave pair- 
ing) [S5]. In this regard, the fact that no strong correlation 
has been observed between the residual resistivity (which in- 
directly characterizes the impurity scattering) and the criti- 
cal temperature of the (nearly) optimally electron-doped 
BKFA (see Table SI) indicates that the level of interband im- 
purities in the Born limit is very small. Thus one only needs 
to estimate the intraband scattering rates 7a, 7b- 



13 meV and 



0.2 -1.7 -1.7 \ 

0.2 -0.25 -0.25 

-5.34 -0.89 0.2 

V -5.34 -0.89 0.2 / 



(SI) 



Negative elements correspond to interband hole-electron re- 
pulsion, while the positive — to intraband attraction. 

In order to apply this full 4-band model to description of 
the transpoit properties one has to take into account an addi- 
tional set of 4 plasma frequencies and a 4x4 matrix of impu- 
rity scattering rates. The latter are difficult to determine the- 
oretically and thus would have to be treated as free pareme- 
ters of the model. It would render the problem highly over- 
parametrized. On the other hand, in the case of BKFA it is 
known that three larger gaps have approximately the same 
value |Aa| = |Ai| w jAs] w IA4I w 9 meV, while the 
smaller gap is |Ab| = IA2I ~ 3 meV. One can assume this 
restriction exact and introduce it into the theory thus reducing 
the original 4-band model to a more tractable 2-band model 
as explained in the following section. 



Reduction to a two-band model 

In general, the superconducting order parameters are a so- 
lution of a linear system of equations 



(S2) 



The Eliashberg functions Bij{uj) satisfy the symmetry rela- 
tions 



N,B,j{lj) ^ NjBj,{lj) 



(S3) 



Theoretical model 

As a starting point we consider a 4— band model based on 
the band-structure calculations with two hole bands and two 
electron bands crossing the Fermi level that has proven suc- 
cessful in accounting for the thermodynamical properties of 
BKFA [S4]. We use the same input parameters, namely, the 
densities of states Ni = 22 Ry-st"\ N2 = 25 Ry-st"\ 
and A^3 = A'4 = 7 Ry-st^^, the first two having a hole 
while the other two an electron character. The main input, 
the spectral function of the intermediate boson, was taken 
following Ref. S9 in the form of a spin-fluctuation spectrum 
Bij{il) — \ijf{Q,/Vlsf ) with a linear cj dependence at low 
frequencies. Here Ay is the coupling constant pairing band i 
with band j and ^sf is a characteristic spin-fluctuation fre- 
quency, the values of which correspond to those in Ref. 4: 



and, therefore, can be represented in the form Bij{uj) ~ 
Uij{Ld)Nj, where Uij is a symmetrical matrix. Further, we 
can construct a functional 

4 4 

5{e J = ^J^? - J2 N.e^U^jN^e, . (S4) 

Equation (S2) then results from minimization of with re- 
spect to Ci. As mentioned above, BKFA has three gaps with 
very close absolute values (the first hole gap has the oppo- 
site sign with respect to the other two). Minimizing func- 
tional Eq. (S4) subject to the additional constraints 63 — 64 — 
— ei — Aa, and 62 = — Ab one finds 
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where the matrix elements satisfy 

iVi (All - 2Ai3 - 2A14) + 7V3A33+iV4A44 



Aab = 



A^i + 7V3 + Na 
N2 (A23 + A24) 



^BA = A23 + A24, 
^BB = A22- 

Assuming the matrix elements Eq. (SI), the following cou- 
pling constants of the reduced 2-band model are obtained: 



^"=1 ^0.5 of ]^I^J = {^^B}. (S6) 



The partial densities of states on the Fermi level of effective 
band A and band B are 

Na = A/"! + A^3 + 7V4 = 36Ry-st"\ 

A^B = A^2 - 25Ry-st-^ (S7) 

Interestingly, even though the main interactions in the 4-band 
model with the coupling constants Eq. (SI) come from the 
nondiagonal elements, in the reduced 2-band counterpart they 
are incorporated into the effective intraband Aaa matrix el- 
ement. Figure S2 presents an interactive simulation of this 
effective two-band model (click on the figure to run the sim- 
ulation and adjust the sliders to set the physical parameters of 
the system). 
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FIG. S2: (INTERACTIVE: click on the image to run the simula- 
tion) Real part of the optical conductivity (top left) and the optical 
scattering rate (top right) of Bao.68Ko.32Fe2As2 in the supercon- 
ducting state at 10 K (blue lines) and the normal state at 40 K (red) 
as obtained experimentally (heavy) and from the effective two-band 
Eliashberg theory (thin), (bottom left) The spin-fluctuation spectrum 
used in the simulation in both the superconducting and the normal 
state, (bottom right) Interactive controls: adjust the sliders to select 
the physical parameters of the model. 
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FIG. S3: Temperature depencence of the (a) free energy and (b) su- 
perconducting gaps the 4-band (lines) and reduced 2-band (symbol) 
models, with coupling matrices Eqs. (S1,S6), respectively. 



Verification of the 2x2 model 

The effective 2-band model closely reproduces all the pre- 
dictions of the 4-band model such as the superconducting 
transition temperature T^. = 38.4 K, superconducting gaps 
Aa = 9.7 meV and Ab — 3.7 meV, free energy and su- 
perconducting gaps as functions of temperature, as shown in 
Figs. S3(a) and (b), respectively. The calculated densities of 
states TVa and A^b are very similar, in accordance with the par- 
tial Sommerfeld constants obtained in the treatment of the spe- 
cific heat data in a phenomenological two-band a-model [S4]: 
7A — 7b (the Sommerfeld constant 7 is related to the density 
of states via 7 = ^iV). 
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